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13.  ABM  1 1 A  C  T 


Sequential  search  procedures  arc  described  for  determining  an 
optimal  dosage  in  the  following  biomedical  problem.  People  are 
subjected  to  a  certain  chcmotheraputic  treatment  and  on  the  one 
hand  It  is  desirable  to  give  each  individual  the  maximal  possible 
dosage.  On  the  other  hand,  high  doses  create  undesirable  toxicity, 
and  it  is  undesirable  to  cross  a  certain  limit  of  allowable 
toxicity  level.  The  optimal  dosage  is  defined  as  the  maximal 
dose  for  which  the  proportion  of  patients  in  the  population  whose 
toxicity  level  will  not  cross  the  allowable  limit  is  y  *  Ih  ~hc 
present  paper  we  discuss  Bayesian  and  non-Bnyesian  sequential 
procedures  for  the  optimal  dosage,  assuming  a  linear  regression 
between  toxicity  and  dosage,  and  normal  conditional  distribution 
of  the  toxicity  level  at  each  dose,  with  a  known  variance.  In 
the  two  models  under  consideration  we  assume  n 

(i)  the  variance  at  dose  x  is  proportional  to  (x~x0)r  for 
x  >  xQ  and  is  zero  for  x  <  xQ;  i' 

(li)  the  variance ■ docs  not  depend  on  the  dose. 
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Sequential  Search  of  an  Optimal  Dosage: 
Some  Preliminary  Results  and  Suggested 
Areas  for  Further  Research,  t) 

B.  H.  Eichhorn  and  S.  Zacks 
Case  Western  Reserve  University 


1.  Introduction. 

The  present  paper  summarizes  our  preliminary  results  in  the  area 
of  sequential  search  for  an  optimal  dosage,  and  brings  forth  several  sug¬ 
gestions  for  further  research.  The  search  procedures  suggested  here  are 
based  on  two  statistical  models  which  are  specified  in  the  next  section. 
These  models  seem  plausible  to  us.  However,  it  may  be  of  great  interest 
to  investigate  how  sensitive  are  the  procedures  to  the  assumptions  of 
these  models.  Suggestions  of  how  to  proceed  in  this  study  of  robustness 
will  be  given  in  Section  7.  Similarly,  wo  may  try  to  change  the  objective 
function  and  investigate  the  possible  implications.  In  Section  2  wc 
specify  the  models  under  consideration  and  the  objective  func Lions. 

Section  3  provides  a  sequential  search  procedure  for  one  of  these  models, 
in  which  the  response  variance  is  proportional  to  the  dosage  squared. 
Section  A  gives  ■  sequential  procedure  for  a  model  of  fixed  response 
variance.  In  Section  0  we  present  Bayesian  sequential  procedures  for  the 
models  mentioned  above.  Monte  Carlo  demonstration  is  exhibited  in  Section 
6.  Section  7  is  devoted  to  open  problems  and  suggestions. 


t)  Partially  supported  by  Project  NR  042-276,  of  the  Office  of  Naval 
Research  at  Case  Western  Reserve  University. 
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2.  The  Statistical  Models. 

Let  x  designate  a  dosage  and  Y(x)  the  toxicity  level  at  x. 
Y(x)  is  a  random  variable.  We  consider  here  the  following  regression 
model.  The  toxicity  level  is  negligible  for  all  0  <  x  <  x^.  From 
the  expected  value  of  Y(x)  is  linear  in  (x-Xq)  with  a  slope  b,  as 
illustrated  in  Figure  1. 

Y 

E{Y(x)}  =  b(x-xQ) 


l 

• 

i 

i 

This  assumption  concerning  the  relationship  between  the  expected  toxicity 
level  e[y(x)}  and  the  dosage  x  is  common  to  all  the  models  specified 
below.  The  assumptions  concerning  the  conditional  distributions  of  Y(x), 
given  x,  and  how  many  of  its  parameters  are  known  constitute  the  following 
two  model ss 

Model  1;  The  conditional  distribution  of  Y(x)  given  x  is  normal  with 

2  2  2 

•  mean  bCx-x^)  and  variance  d  (x-x^)  >  where  o  is  known, 

and  b  unknown.  0  <  b  <  w. 

Model  2:  The  conditional  distribution  of  Y(x)  given  x  is  normal  with 

2 

moan  bCx-x^)  and  a  known  variance,  d  j  b  is  unknown. 
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As  seen  here*  the  present  paper  is  based  on  the  assumption  that  the  value 
2 

of  d  is  known.  Procedures  of  sequential  search  which  are  not  based  on 
2 

a  given  value  of  d  will  be  the  subject  oi  further  investigations. 

Let  T]  designate  the  threshold  of  dangerous  toxicity  levels.  It 
is  desirable  that  Y(x)  will  not  exceed  T}.  Since  Y(x)  is  a  random 

variable  we  cannot  guarantee  with  certainty  that  Y(x)  <  T|  (unless  x  <  Xq 

\ 

We  therefore  require  that  for  some  tolerance  probability  Y»  x  will  be 
such  that  p{Y(x)  <  T]|x}  >  y.  The  largest  value  of  x  for  which  this 
probability  is  at  least  y  will  be  designated  by  5^*  will  be  called 

the  optimal  dose.  This  is  the  unknown  parameter  which  we  wish  to  find  by 
sequential  experimentation.  It  is  easy  to  verify  that  ^  is  given  by 


xQ  +  l/(b+Z  c)  ’  under 


Mrdel  1 


(2.1) 


xq  +  (^"Z^dJ/b  ,  under  Model  2, 


where  Z y  -  &  (y)  1 e  the  y- Tractile  of  the  standard  normal  distribution. 

The  objective  is  to  determine  the  dosages,  x,  at  each  stage  of  experi¬ 
mentation  as  close  as  possible  to  §  and  not  to  exceed  £  .  Thus,  the 
statistical  problem  is  how  to  utilize  the  available  information  on  the  un¬ 
known  parameters  and  determine  a  sequence  of  dosages  x^,  x^,  .  •«,  so  that: 

(i)  For  each  n  =  1,  2,  ... 


(2.2)-  P  fx  <  5  }  >  1  -  a  ,  for  all  0  <  b,  a  <  oo  . 

b,dl  n  -  y  ~  - 

Any  procedure  of  determining  txn$n  >  l}  which  guarantees  (2.2)  is  called 
feasible . 
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(ii)  We  would  further  like  to  guarantee  that  lim  x  =  i  in  prob- 

n.-»  «  n  y 

ability.  Such  a  procedure  will  be  called  const  stent. 

I 

In  order  to  secure  feasibility  and  avoid  certain  theoretical  diffi¬ 
culties  we  assume  a  knowledge  of  a  safj  dose  x*  such  that  x^ < x* < . 
We  can  therefore  restrict  attention  to  procedures  which  assign  dosages  not 
smaller  than  x*  . 


3.  A  Search  Procedure  for  Model  I. 

Given  the  values  of  x^,  • 

Y,»  ...»  Y  we  determine 
1*  '  n 


(3.1) 


i  n 

=  i  E 

0  i=l 


x  .  and  the  associated 
n’ 


Y  values 


Let  U+  =  max(0,Un).  The  value  of  xn+1  is  determined  then  as  a  function 
of  Un  according  to  the  formula: 


(3.2) 


We  prove  now  that  this  procedure  is  (i)  feasible;  (ii)  consistent  and 
(iii)  optimal  in  a  certain  sense. 


(i)  First,,  consider  the  distribution  of  =  Y^/(x^-xQ),  i  =  1,  ..., 

Since  x^  is  a  function  only  of  (Y^,  ...,  Y^  ^),  the  conditional  distribu 

tion  of  U.  given  the  c-field  (Y^,  •  ••>  Y^_,)  is  normal  with 

2  2 

expectation  b  and  variance  c  .  Since  b  and  a  are  independent  of 

we  obtain  that  is  independent  of  ^  and  has  a  normal  distri- 

2  ’  2 
bution  7({b ,a  ).  Hence,  U^,  ...,  are  i.i.d.  with  ?i[(b,d  )  distribu¬ 


tion,  This  implies  that 
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(3.3)  p{*0  t  „/(i£  +  a  (-^  +  Zy))  <  Sy  }  = 


p{uh  +  —  z  >  0}  >  p{u  +  —  z  >  b) 
ln  JZ  1"“  -  ;  '  ;  n  "  ; 

Let  *„  =  {x0  t  ,/(ft  +  a(^?  +  Zy))  <  *g}  . 

Q 

and  let  A  designate  the  compliment  of  A  .  Then, 
n  n 


=  1  - 


*n+l  =  x*  ?W  +  [o  +  l/(^  +  a(-j^  +  Zy))]  1  tAn) 
where  l{*}  is  the  indicator  function.  Finally,  since  x*  <  |  , 

Pk-H  &  5y}  2  p{><0  +  ’  A<  +  «£jf  +  Zy))  2  5V}  • 

This  proves  that  the  procedure  is  feasible. 


(ii)  By  the  strong  law  of  large  numbers,  -♦  b  with  probability  one. 

Hence,  since  b  >  0  xn  -♦  %  with  probability  one.  This  proves  the  (strong) 
consistency  of  the  procedure. 

(iii)  The  procedure  prescribed  by  (3.1)  and  (3.2)  has  the  following 
optimality  property: 

Among  all  feasible  procedures  the  procedure  [xn!n  >  l}  proscribed 
bx  (3.1)  -  (3.2)  is.  asymptotically  uniformly  most  accurate,  i.e.,  if 
[xn5n  >  l]  is  any  sequence  of  dosages  which  satisfies  the  feasibility 
condition  (2.2)  then 

PbJ'*n+l  <  *  Pb,oK+l  <  «')  ** 

all  (b,a)j  and  every  §’  <  5^  »  a ni^  n  sufficiently  large 
(may  depend  on  5’)* 
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The  proof  of  (3.4)  is  based  on  the  fact  that  b  =  U  +  —  Z, 

c,n  r»  ^  1-a 

is  a  uniformly  most  accurate  upper  confidence  limit  for  b  (see  Lehmann  * 

A 

pp.  78-81)  for  each  n  =  1,  2,  ...  .  Hence,  if  b  is  any  other  upper 

a  >  n 

confidence  limit  for  b  then 


(3.6) 


Pb„A,n  >  b,S  i  Pb,oSSa,n  >  >''!  for  cach  ,  =  1' 
and  all  (b,o)  and  every  b'  >  b. 


Let  b*  =  U*  +  ~  Z,  then  for  n>( — we  have  fb*  >b'}c[b  >b'l.  Letting 

a,n  n  n  1-a  \  b,  •  4  a,n  J  a,n 

b  =  n/(x  ii  "  xr )  -  °  Z  we  obtain  from  (3.b)  and  (3.2)  that  (3.4) 
cl  }  n  v  n  *  l  u  Y 

holds,  for  these  values  of  n. 


4.  A  Search  Procedure  for  Model  7'.. 

Under  Model  II  there  are  several  complications.  We  still  wish  to 
estimate  b  by  U^,  as  given  by  (3.1).  However,  since  the  sequence  of 
values  is  not  determined  before  the  observations  commence,  the  sampling 
distribution  of  is  not  normal.  We  cannot  even  apply  the  usual 

versions  of  the  Central  Limit  Theorem,  since  U^,  U^,  ...,  are  dependent. 
It  is  easy  to  prove  that  E|Un]  ~  b  for  all  (b,o),  and  furthermore 
cov(U.,Uj)  -  0  for  all  1  <  i  <  j  <  n.  Indeed,  for  all  i  <  j, 


(4.1) 


E^U.}  =  e[u.e{u  ./*.}}  = 


Hence  U^  and  are  uncorrclated.  From  this  result  we  immediately 


imply  that 


(4.2) 


Var(ll  }  =4  £  E / - 1 - 

n2  i,l  \(x...x0)2/ 


E. Lehman,  Testing  Statistical  Hypotheses. 
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We  notice  in  this  formula  that  Var{Un]  may  diverge  to  infinity  if 

E{(x^-Xq)  }  -  •  In  order  to  avoid  such  a  possibility  we  require  that 

all  will  be  greater  or  equal  to  x^,  whore  xQ  >  x^.  Then, 

2 
o 

n(xo‘xo} 

Thus  we  consider  the  following  procedure:  After  observing  Y^,  Yn 

1  n 

at  x.,  ....  x  compute  U  -  —  £)  Y,/(x.-x_)  and  set 
1  n  r  n  n  r  i  O' 

(4.4)  xn+1  =  max^tj,x0  +  (T|-?yO/^Un  +  *|  dS^  , 

where  d  is  determined  so  that  x  ,,  is  a  feasible  and  consistent 
a  n+1 

sequence.  Employing  the  Chebytchev's  inequality  one  could  use,  for  example 
d^  -  a  8  *  ^his  is  however  a  very  conservative  approach.  There 

is  no  doubt  that  more  refined  bounds  could  be  determined,  so  that  a  faster 
convergence  could  be  expected. 

5.  Bayes  Procedures  for  Mode  1  I  and  Model  II . 

In  a  Bayesian  framework  we  assume  that  the  unknown  parameter  is  a 
random  variable  having  some  prior  distribution.  After  observing  n  obser¬ 
vations  (n  =  1,  2,  ...)  of  Y^j  ...,  Yn  at  dosages  x^,  ...,  we 
convert  the  prior  distribution  of  b  to  a  posterior  distribution  given 
5"^.  We  then  determine  x  so  that  the  posterior  probability  of 

fx  ,,  <  E  )  given  tf  will  be  at  least  1  -  a.  This  will  assure  that  the 

t  ni  l  -  J  n 


(4.3) 


Var[U  1  £ 
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expectation  of  the  left  hand  side  of  (2.2)  with  respect  to  the  prior 
distribution  of  b  will  be  at  least  1  -  a«  This  property  docs  not 
imply  the  feasibility  condition  (2.2).  The  property  that  we  attain  is 

somewhat  weaker  and  actually  depends  on  the  prior  distribution  assumed. 

We  shall  therefore  say  that  a  procedure  is  Bayes-feasible  with  respect  to 
a  prior  distribution  H  of  b  if 

S 

(5.1)  PH^Xn+l  -  -  1  "  a  for  overy  n  =  2»  ••• 


and  all  d,  0<n<®,  pH^xn+l  -  designates  the  posterior 

probability  with  respect  to  the  distribution  of  b. 

We  provide  now  explicit  formulae  for  the  determination  of 
{xn+i,  n  =  1,  2,  ...}  satisfying  (5.1),  for  a  normal  prior  distribution 
of  b,  ??(P,Vq)j  with  prior  moan  p  and  prior  variance  VQ.  We  remark 
here  tiiat  one  could  use  the  same  methodology  to  derive  proper  formulae 
for  other  prior  distributions  of  b.  Let  X  =  xn  “  xq>  n  ~  ••• 

and  as  before  U.  ~  Y./X.,  i  =  1,  ...,  n.  To  give  a  general  framework 

t 

d^X^  ,  under  Model  1 

n  ’ 

2 

0  ,  under  Model  II  . 


Given  Y^,  ...,  Yn  and  X^,  ...,  Xf),  it  is  easy  to  prove  that,  if  b  has 
a  prior  normal  distribution  ‘//(j3,Vq)  then,  its  posterior  distribution  is 
also  normal,  '/?(pn,Vn),  with  posterior  moan 
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(5.3) 


A  =  11  .  4-  (y  -  fl  .X  )X  V  ./(t *“  +  X2V  ,),  n  -  1,  2, 
*n  *n-l  '  n  *n-l  n'  n  n-1' '  n  n  n-1 *  * 


and  posterior  variance 


(5.4) 


Vn  ^  Vl 


n 


T2  +  X2  V 
n  n  n-1 


n  =  1,  2, 


From  the  recursive  formulae  (5.4)  we  obtain  that  the  posterior  variance 
of  b  is 


Thus,  if  all  It  >  X*  then  is  in  order  of  magnitude  (in  probability) 

of  n  f  i.o. ,  Vn  =  Op(n  ^)  as  n  -♦  »  . 

Similarly  we  obtain  that  the  posterior  mean  pn  is  given  explicitly, 
under  Model  Iv  as: 


(5.6)  pn  =  Un(H  d2/nV0)_1  +  pQ(l +o2/V0)/(n  td2A0)  . 

Hence,  under  Model  I,  *♦  b  with  probability  one,  as  n  -♦  «  .  An  explicit 
formula  of  p^  for  Model  II  is  considerably  more  complicated  and  we  shall 
use  the  recursive  formula  (5.3)  with  =  a  .  It  can  be  shown  that  p 
is  a  consistent  estimator  of  b,  for  almost  every  b  (with  regard  to  H) 


alsci  under  Model  II. 
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The  Beyer,  procedure  specifier,  the  following  sequence  of  dosages: 


xn+l  s  maxK»6n,Y^  WhcrG 


(b.  7) 


xO  +  1l/tpn  +  +  ^1_a  ^  ]  ,  for  Model  1 

XQ  1  (11  -  2Yd)/(pu  +  Z:_q  ,/T)  ,  for  Model  11  . 


6.  Monte  Carlo  Comparisons. 

In  the  present  section  we  compare  the  vaiious  sequential  procedures 
numerically  by  starting  with  an  initial  dose  x^  simulating  deter¬ 

mining  *2  simulating  Y2  etc.  We  present  the  results  of  50  such  itera¬ 
tions,  The  parameters  of  this  simulation  a:o: 

b  =  3. ,  Xq  —  O.j  Xq  —  1  •  j  T|  -  10. ,  0  ■*  l,j  a  =  .05$  Y  “  .99*- 

The  initial  dosage  Is  x^  -  3,5.  The  value  of  £  i  s 

^  1.878  ,  for  Model  1 


l  2.558  ,  for  Model  II  . 


V 


0 


.25 


p  —  2. 86 


Table  1.  Simulated  Dosage  Determination 
and  Baycj,  (lj.)  Procedure!. 


for  tile  non-Bayun  (N.H. ) 


_  Model  1 

n  N.n.  b. 


Mod  e J  ft 

N.  B.  B. 


1 

2 
3 
6 
8 
6 
7 
is 
9 
10 
I  1 
12 
13 
16 

15 

16 

17 

18 

1  9 
20 

2  J 
22 
23 

2/i 

20 

26 

27 

28 

29 

30 

31 

32 

33 
3/1 
30 

36 

37 

38 

39 
60 
6  J 
6ft 
63 
66 
68 
66 
67 
6  8 
6  9 
SO 


1.166  1 

1 • 8678 

1  .2078 

1 . 870 J 

1 • 3668 

1 . 8737 

1 . 6273 

1 . 8903 

1 • 6266 

1 . 87  78 

1 .6983 

1.61 22 

1  •  3  3  7  3 

i  . 6328 

J  .  88  89 

1 . 6398 

I  .  62  38 

1 . 6867 

1 .6196 

1 . 6779 

1  •  66  27 

1 . 6928 

1  •  68  83 

1 . 6  8  6  3 

1 .6609 

1 . 6863 

1 . 66  60 

1  *6807 

1 . 6986 

1 • 7296 

1 .6966 

1  • 7880 

1 .6896 

1.71 99 

1 . 7081 

1.781 8 

1 . 7367 

1 • 7860 

1  •  76  86 

1  • 766  6 

1.7221 

1 . 7638 

1.7181 

1 • 7397 

1.7119 

1 . 7336 

1.71  19 

1 . 7330 

1  •  7298' 

1 .768  1 

J . 7822 

1  •  7676 

1.7613 

1 • 7  78  1 

1.77 8  8 

1 . 7870 

1  •  7790 

l  •  7902 

1.776ft 

1 . 7  68  7 

1  •  76  66 

1  •  7  7  86 

1  •  7788 

1 .7866 

I  .  7766 

1  •  7  882 

1.7712 

1  .  7  820 

1 . 7  67  8 

1.77 1  7 

1 . 76  1  8 

I  • 7  730 

I • 7637 

1 . 7  76  3 

1.73  9  7 

1 . 7  7  07 

J . 7660 

1 . 7  7  66 

1.7639 

1.7761 

1 . 7609 

1.7712 

1 . 7767 

1 . 7 887 

1 . 7636 

1  •  7  9  I  7 

1.7917 

1 . 799 1 

1  •  8  0  2  6 

1 . 8093 

1 • 3096 

1 . 6183 

1.61  /;  1 

1 . 6190 

1  •  607  6 

1 . 8 1 80 

1 . 7966 

1 . 6066 

1 • 6063 

1.61  60 

1 . 000  0 

ft.  0660 

1  •  1  066 

ft.  03  13 

1  • 22 1  1 

ft.  0763 

1.3313 

Pr  1  ft!!  1 

1  •  3683 

ft  .  1  J  0  1 

1  • 6  366 

2. 1682 

1 . 8206 

ft. 2010 

1 . 86  1  7 

ft  •  2  1  7  1 

1  • 66  20 

ft. 2827 

1 .60  00 

ft.  2787 

1 *6986 

ft. 3027 

1.7) 0 J 

ft. 2983 

1  • 7  303 

ft. 30 12 

1.7673 

ft. 3067 

’ • 8068 

ft. 36  18 

1 . 8166 

ft. 3380 

1 . 6266 

ft  *  3836 

1 . 8686 

ft. 3702 

1  *  86  36 

ft. 6021 

1 • 9027 

ft.  6 1 62 

1 . 8928 

2 . 3693 

1 .8983 

2. 3886 

J  .  9012 

2. 3790 

J . 9090 

ft. 3793 

1 .9316 

2. 3986 

1  •  9876 

2.6 227 

1  •  9  7  22 

2.6328 

1.9909  ' 

ft.  66  88 

1 . 99  99 

ft.  6621 

2.  002 1 

2.  66  72 

ft.  00  ]  6 

2.639ft 

ft • 0 1 30 

2. 66  93 

ft. 01 92 

2. 6679 

ft.  08  1  7 

2. 66  6  7 

ft.  0238 

2. 66  1  2 

2. 0236 

2. 63 0A 

ft • 0297 

ft  .  6  3  7  1 

ft. 0307 

2. 633 1 

ft. 0387 

2.6378 

2*  0  6  2  6 

2. 6377 

ft*  066  1 

2. 6267 

ft. 0698 

2.  669  3 

ft. 07  1  2 

2. 6  SOB 

ft.  06  1  0 

2. 66 78 

ft ■ 096/ 

2.67  ) 

ft  .  1  Ci  2  9 

2.  6 Si  6 

ft .  11  1  3 

2 .  6  9  1  6 

2. 1 096 

2.6  86  3 

ft  .  1  f-  j  7 

2. 6730 

ft  .  ]  1  0  3 

2.  /.; 

-12- 


The  simulated  dosage  values  which  arc  exhibited  in  Table  1  illustrate 
the  approach  of  these  sequences  of  doses  to  the  optimal  doses,  given  by 
5^.  In  Model  I  the  Bayes  procedure  yields  dosages  which  arc  somewhat 
closer  to  the  optimal.  The  differences  between  the  non-Bayes  and  the 
Bayes  procedures  become  insignificant  as  n  grows.  This  phenomenon  depends 
however  on  a  "good"  choice  of  prior  parameters  for  the  Bayes  procedure.  A 
similar  phenomenon  is  observed  in  Model  II.  The  non-Bayes  procedure  (de¬ 
rived  in  Section  4)  yields  at  the  beginning  values  of  x  close  to  x* 
and  the  convergence  is  slow.  This  is  due  to  the  "over  pessimistic"  choice 
of  the  parameter  in  (4.4).  We  further  observe  that  in  no  case  the 

dosages  obtained  exceed  the  optimal  dose  This  is  valuable  character¬ 

istic  of  the  proposed  procedures. 

7.  Suggestions  for  Further  Hcsoarc.h. 

The  following  suggestions  for  further  research  arc  based  on  our  own 
speculations  concerning  the  relevancy  of  the  models  studied  in  the  previous 
sections.  Our  list  of  open  problems  is  classified  into  two  major  classes, 
(i)  Variations  of  the  statistical  assumptions  and  objectives  within  the 
framework  of  the  present  problem,  (ii)  Extension  of  the  models  info 
multivariate,  multi -dimensional  and  time  dependent  problems.  Within  the 
first  major  class  wc  suggest  to  consider  the  following  problems 
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(1.1)  The  conditional  distribution  of  toxicity  levels,  Y(x),  for  a 
given  dosage,  x,  is  normal  with  mean  bCx-x^)  and  unknown  variance 
a  (x).  This  variant  of  the  statistical  model  can  be  further  broken  up 
into  several  special  cases. 

(1.2)  Problems  connected  with  unknown  intercepts  (place  of  Xq)  and 
known  or  unknown  variances,  assuming  still  normal  conditional  distribu¬ 
tions  . 

(1.3)  The  effect  on  the  procedures  caused  by  deviations  from  normality  of 
the  conditional  distributions  around  the  regression  lines. 

(1.4)  derivation  of  search  procedures  when  the  toxicity  can  assume  values 
only  on  a  discrete  sot. 

(1.5)  Deviations  from  linearity  of  the  toxicity-dosage  regression  line. 

(1.6)  Sensitivity  analysis  -  study  of  the  robustness  of  the  search  pro¬ 
cedures  concerning  the  basic-  assumptions  on  the  distributions  and  on  the 
toxicity-dosagc  regression. 

(1.7)  Formulation  of  different  types  of  objective  functions. 

In  the  second  class  of  open  problems  wc  mention: 

(2.1)  Multivariate  response  -  the  observations  consist  of  vectors  of 
several  components,  one  of  which  is  toxicity. 
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(2.2)  Time  dependent  problems  -  patients  are  subject  to  continual  treat¬ 
ment.  The  effect  of  prior  treatments  on  future  dosages,  employing  the 
information  gathered  on  each  individual  separately.  The  determination 

of  the  optimal  spacing  between  epochs  of  treatments, 

(2.3)  The  search  for  the  optimal  combination  of  various  drugs. 
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